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Abstract. Using the IRAM 30 m telescope, we report the detection of the CO(3-2), CO(4-3), CO(5-4) and CO(6- 
5) lines in the gravitational lensed submm galaxy SMM J16359+6612 at 2; = 2.5. The CO lines have a double 
peak profile in all transitions. From a Gaussian decomposition of the spectra we show that the CO line ratios, 
and therefore the underlying physical conditions of the gas, are similar for the blue and the redshifted component. 
The CO line Spectral Energy Distribution (SED; i.e. flux density vs. rotational quantum number) turns over 
already at the CO 5-4 transition which shows that the molecular gas is less excited than in nearby starburst 
galaxies and high-z QSOs. This difference mainly arises from a lower average H2 density, which indicates that 
the gas is less centrally concentrated than in nuclear starburst regions in local galaxies. We suggest that the bulk 
of the molecular gas in SMM J16359+6612 may arise from an overlap region of two merging galaxies. The low 
gas density and clear velocity separation may reflect an evolutionary stage of the merger event that is in between 
those seen in the Antennae and in the more evolved ultraluminous infrared galaxies (ULIRGs) like e.g. Mrk231. 
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1. Introduction 

The intensity of the far-infrared (FIR) background 
indicates that emission from warm dust contributes 
significantly to galaxies' overall energy output over the 

mm 



Hubble time (Puget etal. HMHl Fixsen etal. T^). 
Recent surveys at submm and mm wavelengths with 
SCUBA at the JCMT and MAMBO at the IRAM 30 m 
telescope have identified part of the underlying galaxy 
population responsible for the strong FIR emission. Up 
to now several hundred submm/mm galaxies (SMGs) 
have been identified (e.g. Small et al. 19T)T, Bertoldi et al. 
IMTOl Ivison etal. Webb etal. 2003, Greve etal. 

I2004|l . These surveys, however, are limited by confusion 
and therefore only trace the bright part of the SMG 
luminosity function (Blain etal. 2002). Observations of 
CO, providing valuable information on the dynamics, 
size and mass of the molecular reservoirs in these objects, 
have only been reported for 12 SMGs so far (e.g. Frayer 
et al. [TUMI Ivison et al. I^UHTl Downes & Solomon I^DUSl 
Genzel et al. 11003 Neri et al.EOOSI Greve et al. OnnSjl . Due 
to their selection based on the mm/submm continuum, 
all these sources, however, are intrinsically luminous, and 
not representative of the faint end of the SMG luminosity 
function that dominates the FIR background. 



Recently Kneib et al. H2004|l discovered the strongly 
lensed submillimeter galaxy SMM J16359+6612 towards 
the galaxy cluster A 2218. This source is lensed by the 
cluster into 3 discrete images and the large magnification 
(14, 22 & 9 for images A, B & C respectively) implies 
that its intrinsic submm flux density (S'sso/^m = 0.8 mJy, 
corresponding to Lfir « 10^^ (Kneib et al. l2UUHl ') 
is below the confusion limit of existing 850 ^m surveys 
(« 2 mJy, Blain et al. l2002|l . It therefore provides a unique 
opportunity to investigate a source which is presumably 
more representative for the submm population. 
Detections of CO towards SMM J16359+6612 have been 
reported by Sheth et al. I|2004l CO 3-2) and Kneib et al. 
proSl CO 3-2 & 7-6). In this letter we report on the 
detection of the CO 3-2, 4-3, 5-4 & 6-5 lines towards the 
strongest lensed component B. For the derived quantities 
in this paper, we use a A cosmology with Hq — 71 
kms"iMpc~\ = 0.73 and n,n = 0.27 (Spergel etal. 

mm. 



2. Observations 

Observations towards SMM J16359+6612 B were made 
with the IRAM 30 m telescope during nine runs between 
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Fig. 1. Spectra of the CO(3-2), CO(4-3), CO(5-4) and CO(6-5) lines towards SMM J16359+6612 B, with Gaussian 
fit profiles superposed. The velocity scale is relative to a CO redshift of z = 2.5174. The velocity resolution is 50 kms"^ 
(3-2, 4-3, 5-4) and 55 kms~^(6-5). All spectra are shown with the same flux scale. 



Jan. and March 2005 in mostly excellent weather con- 
ditions. We used the AB and CD receiver configuration 
with the A/B receivers tuned to the CO(3-2) (98.310 GHz, 
3 mm band) and CO(6-5) (196.586 GHz 1mm band) and 
C/D to CO(4-3) or CO(5-4) (131.074, 161.637GHz, 2 mm 
band). The beam sizes/antenna gains for increasing ob- 
serving frequencies are 25"/6.1JyK~^, 19"/6.5 JyK~^, 
15'76.9JyK'i and 12.5"/7.7 Jy K^i. Typical system 
temperatures were «120K, 220 K, and 270 K (T^) for 
the 3, 2, and 1 mm observations. The observations were 
done in wobbler switching mode, with switching frequen- 
cies of 0.5 Hz and a wobbler throw of 50" in azimuth. 
Pointing was checked frequently and was found to be sta- 
ble to within 3". Calibration was done every 12min using 
the standard hot/cold-load absorber measurements. The 
planets Uranus and Neptune were used for absolute flux 
calibration. We estimate the flux density scale to be accu- 
rate to about ±10-15%. 

Data were taken with the IMHz filter banks on the A/B 
3mm receivers (512 channels, 512 MHz bandwidth, IMHz 
channel spacing) and the 4 MHz filter banks for the 2 and 
1.3 mm observations (256 channels, 1 GHz bandwidth, 
4 MHz channel spacing). The data were processed with 
the CLASS software. We omitted scans with distorted 
baselines, and only subtracted linear baselines from in- 
dividual spectra. Their average was then regridded to 
a velocity resolution of 50 kms~^ (CO 3-2, 4-3, 5-4 ) 
and 55 km s"^ (6-5) leading to rms noise values (T^) of 
0.26mK (1.6mJy), 0.47mK (3.0mJy), 0.6mK (4.1mJy) 
and 0.38 mK (2.9 mJy) respectively. The total observing 
time in the final spectra is 8.7h, 6.8h, 4.7h and 8.2h for 
the 3-2 to 6-5 lines, respectively. The final spectra are 
presented in Fig. We also generated an averaged CO 
spectrum by averaging all individual CO transitions, with 
equal weight. This spectrum is shown at a velocity reso- 
lution of 15 kms^^ in Fig.[21 

3. Results 

We detected all observed CO transitions from 
SMMJ16359-h6612 B. The hne profiles for all four 
lines are similar, and show the characteristic double peak 
(see Fig.^ already recognized in previous interferometric 
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Fig. 2. Averaged spectrum of all four CO lines at a veloc- 
ity resolution of 15 kms~^ with Gaussian fit profile super- 
posed. All CO lines have equal weight. The velocity scale 
is relative to z = 2.5174. The total integration time for 
the spectrum is 28.4 hours (on + off) 

studies (Sheth et al. I^UUl Kneib etal. lM75|l . 

The small separation between the lensed images A and 
B of 15.0" (Kneib etal. EEHI) implies that our spectra 
also have a contribution from component A. Taking the 
relative integrated flux densities between component 
A/B and the beam size of the 30 m telescope at the 
four observing frequencies into account, we estimate the 
contribution from component A to the integrated line 
intensities of our spectra to be 24%, 12%, 4% and 2% for 
the 3-2 to 6-5 transitions respectively. These corrections 
have been applied to the following quantitative analysis. 
Due to the larger distance to component C and its lower 
magnification, the contribution from this component is 
negligible. 
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Table 2. CO line parameters towards SMMJ16359+6612 B, corrected for the contribution from image A. The mean 
redshift is the flux-weighted average of the redshift of the red (i?) and blue (B) component. SVr^b denotes the peak- 
to-peak velocity separation between the red and blue components. Parameters have been derived from Gaussian fits 
by keeping the line widths fixed at 250 and 160 kms~^ for the blue and red components. 



Transition 


■2' mean 


S„n 


Si, B 




I 


Ib/Ir 


L' " 








[mJy] 


[mJy] 


[kms~^] 


[Jy kms""""] 




[10^ K kms"^ pc^] 




CO(J=3^2) 


2.51737 


8.5 ±0.7 


5.3 ±0.6 


310 ± 20 


2.8 ±0.2 


1.0 ±0.1 


4.4 ±0.3 


5.8 ±0.4 


CO(J=4^3) 


2.51732 


11.9 ± 1.5 


7.3 ± 1.2 


340 ± 25 


4.0 ± 0.4 


1.0 ±0.2 


3.4 ±0.3 


10.7 ± 1.1 


CO(J=5^4) 


2.51746 


15.1 ± 2.0 


9.7 ± 1.6 


295 ± 30 


5.1 ±0.5 


1.0 ±0.2 


2.8 ±0.3 


17.3 ±1.8 


CO(J=6^5) 


2.51738 


12.5 ± 1.9 


7.0 ± 1.5 


340 ± 35 


4.0 ± 0.5 


0.9 ±0.2 


1.5 ±0.2 


16.1 ± 2.1 



" corrected for a magnification m — 22 (Kneib et al. l2UU4t 
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Fig. 3. Left: Observed, velocity-integrated CO fluxes vs. rotational quantum number (CO hne SED). Filled squares 
are the integrated flux densities listed in Table[21 open squares were adopted from Kneib et al. I|2005|l . The upper three 
curves are LVG models for the entire line profile (red+blue components added together). The lines show the LVG- 
predicted fluxes for ^(Hz), Tki„ of lO^ ^ cm'^, 35 K (solid), lO^ i cm'^, 80 K (dashed) and lO^-^ cm'^, 20 K (dotted). 
The lower data points are the fluxes for the red (open triangles) and blue (open circles) components taken separately. 
The lower curve shows the same LVG model as the solid upper curve. Right: Comparison of various normalized CO 
line SEDs: SEDs are shown for BR1202-0725 {z = 4.7, open circles Carilli et al. l20()2jl . the high-excitation component 
in the center of M82 (WeiB et al. I2005L no data points plotted), NGC 253 center (filled squares taken from Bayet et al. 
I2004|l . SMM J16359+6612 (no data points plotted), the overlap region in the Antennae galaxies (open triangles Zhu 
et al. l2()()3|l and the inner disk of the Milky Way (filled triangles, Fixsen et al. 1199911 . The CO line SEDs are normalized 
by their CO(l-O) flux density. 



Since the line profiles for all transitions are similar, we 
used the average spectrum of all CO lines (Fig.jSJ to 
determine the line parameters of the blue and redshifted 
component in the spectrum. A Gaussian decomposition 
yields line widths of 250 ± 45kms-^ and 160 ± 20kms-^ 
for the blue and the red components respectively. The 
integrated flux density ratio between the blue and the 
red component is close to unity {Ib/Ir = 0.9 ± 0.2). The 
parameters of the spectral decomposition are summarized 



in Table U] 

For the Gaussian decomposition of the individual CO 
lines we fixed the line width in both components to 
the values derived above. The line parameters for each 
transition are given in Table |21 Our integrated CO (3-2) 
line flux is in agreement with the Bure flux reported by 
Kneib etal. H2005fl (2.8 vs 2.5 Jy kms"^) but lower than 
the flux reported by Sheth etal. (3.5Jy kms-\ I^UUljl . 
The integrated and peak flux densities of the CO lines are 
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Table 1. Gaussian decomposition of the averaged CO line 
spectrum (see FigEJ- 





blue 


red 


z 


2.5155 ±0.0002 


2.5192 ±0.0001 


5v [kms"^]" 


-161 ± 17 


±156 ±9 


FWHM [kms~^] 


250 ± 45 


160 ± 20 


Su norm 


0.58 ±0.18 


1 


-^norm 


0.91 ±0.16 


1 



Center velocity relative to z=2.5174 (Kneib etal. l2U05t . 



rising with increasing rotational quantum number up to 
the 5-4 line. The line luminosities (L', see e.g. Solomon 
et al. H1997|l for the definition) show that already the 4-3 
line is subthermally excited relative to the 3-2 line. 
To investigate the CO excitation in more detail we used 
a spherical, one-component, large velocity gradient model 
(LVG). We use the collision rates from Flower H2(J(J1|I with 
an ortho/para H2 ratio of 3 and a CO abundance per 
velocity gradient of [CO]/dv/dr = 8 • lO'^ pc (kms"^)-^ 
All lines are well fitted with a single LVG model. A 
good fit to the observations is provided by n(H2) = 
lO^-'^ cm~'^ and Tkin = 35 K. Other temperature-density 
combinations with similar H2 pressure also match the 
data (see Fig.lSJ. H2 densities above n(H2)= lO"^-^ cm"-^, 
however, do not provide a good fit to the observations. 
In Fig.Ol we compare selected LVG models with the 
observed integrated flux densities including the CO 
3-2 and 7-6 measurements by Kneib etal. 1)2005(1 . 
We estimated the 7-6 flux using a line ratio of 
^co(7~6)/-^co(3-2) ~ 0.7 (line integral corrected for 
the continuum emission at 1mm). The flgure shows that 
the predicted integrated flux density for the CO (1-0) 
line is not a strong function of the model parameters. 
Our models predict /co(i-o) ~ 0.3 — 0.4 Jy kms^^ or 
-^CO(i-o) ^ ^•-'^ ~ ' 10^ K kms^^ pc^ (corrected for a 
magnification m = 22). 

Our Gaussian decomposition of the individual CO transi- 
tions shows that the integrated flux density ratio between 
the blue and the red part of the spectrum is near unity 
for all CO transitions (see Ib/Ir in Table|21l. This implies 
that gas masses, CO line ratios and also the underlying 
excitation conditions are similar in the red and the blue 
components. 

4. Discussion 

At first glance the double horn line profile and the sim- 
ilarity of the gas properties in both spectral components 
suggests that we are looking at a circumnuclear molecular 
toroid. However, if such a stable molecular gas distribution 
would supply the large star formation rate (« 500 Mq yr~^ 
Kneib et al. I2004|l we would expect much higher gas exci- 
tation (the same is true if we were looking at independent 
molecular disks surrounding the nuclei of two merging 



galaxies). This is exemplified in Fig. |2| where we compare 
the normalized CO line SEDs (flux density vs. rotational 
quantum number) of SMM J16359±6612 to those in other 
starburst galaxies/QSOs. The normalization is relative to 
the integrated CO(l-O) flux density (and therefore to the 
molecular gas mass) . These CO SEDs can be used to con- 
strain the CO excitation as the turnover of the CO line 
SED is a sensitive measure of the underlying gas density 
and temperature. 

In SMMJ16359±6612 the CO SED turnover occurs at 
the 5-4 line, which shows that the average gas excita- 
tion is lower than in the center of the well-studied nearby 
starburst galaxies NGC253 and M82. In both cases, the 
CO SED turns over at the 6-5 fine (Bradford et al. IMM 
Bayet etal. 120031 Mao etal. lMTUl WeiB et al. "2005 ) . Even 
higher CO excitation is present in Henize2-10 (Bayet 
et al. IMlljl and some high-z QSOs such as BR1202-0725 
(z = 4.7, Omont et al. ll^Ml Carilh et al. HUU^ and the 
Cloverleaf (z = 2.6, Barvainis et al. I19i?7|l where the flux 
density is still rising even up to the CO(7-6) transition. 
This comparison shows that the molecular gas pressure 
in SMM J16359±6612 is lower than in other active sys- 
tems. From our models, we conclude that the difference 
mainly arises due to the moderate average gas density 
(n(H2) < lO^cm-S) in SMM J16359±6612. This suggests 
that the molecular gas is less centrally concentrated than 
gas in starburst nuclei and QSO host galaxies. This view 
is supported by the large extent (1.5-3 kpc) of the molec- 
ular emission region found by Kneib etal. ((2003 ■ 
One possiblity to explain the combination of lower CO 
excitation and high FIR luminosity {IO^^Lq, Kneib etal. 
12004(1 may be that the CO flux comes from an 'over- 
lap' region between two merging galaxies. Such a geome- 
try has first been suggested by Kneib et al. based on the 
UV/optical morphology of the potential two nuclei (their 
regions a and /?) and the reddened overlap region (their 
region 7). Such active overlap regions are often found in 
nearby interacting systems (e.g., the Antennae: Wilson 
et al.EOOOl NGC 6090: Bryant & Scoville lTM??! Wang et al. 
.20041 . 

For comparison, we therefore also show the CO SED for 
the overlap region in the Antennae in Fig. |3| (Zhu priv. 
com.). Although no observations above the CO(3-2) line 
have been published in literature so far, the low average 
CO(3-2)/CO(l-0) line ratio suggests that the gas in the 
Antennae overlap region has a lower excitation than the 
starburst galaxies/QSOs and SMM J16359±6612. More 
detailed models of the molecular gas in the Antennae show 
that the bulk of the CO emission arises from gas at moder- 
ate density of n(H2)« 10^ cni"^ (Zhu et al. I^UnSjl . Given 
the early stage of the merging process of the Antennae 
and its low FIR luminosity (~ 510^°Lo, Sanders etal. 
nnn^ we suggest that SMM J16359±6612 may be a more 
advanced merger (with higher FIR luminosity and higher 
mean gas density than in the Antennae overlap region) 
but not yet in the typical ULIRG stage (with the gas con- 
centrated in the central < 1 kpc region), such as, e.g., 
Mrk 231. 
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